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(57) Abstract: A method for measuring the temperature rise in anatom- 
14 ical tissue as a result of ultrasound treatment. A first ultrasound signal 

is obtained prior to treating the anatomical tissue, then a second ultra- 
sound signal is obtained after the tissue is treated. Complex analytic sig- 
nals are computed from the first and second ultrasound signals, then the 
depth-dependent delay is computed from the complex analytic signals. 
An echo strain map is generated from the slope of the depth-dependent 
-18 delay. The echo strain map is used to estimate the amount of tempera- 

( ture rise from the first ultrasound signal to the second ultrasound signal. 

An image may then be created showing where temperature rise is occur- 
ring in the anatomical tissue. 
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METHOD FOR MAPPING TEMPERATURE RISE 
USING PULSE-ECHO ULTRASOUND 



Related Applications 

This is a continuation-in-part of Application No. 10/153,241, filed May 22, 
2002, which claims priority to US Provisional Application Serial Number 60/294,135 
filed May 29, 2001. The entire disclosures are incorporated herein by reference. 
Field of the Invention 

[0001] The present invention relates generally to ultrasound, and more particularly, to 
a method for measuring temperature rise using pulse-echo ultrasound waves. 
Background of the Invention 

[0002] Ultrasound medical systems and methods include ultrasound imaging of 
anatomical tissue to identify tissue for medical treatment. Ultrasound may also be used 
to medically treat and destroy unwanted tissue by heating the tissue. Imaging is done 
using low-intensity ultrasound waves, while medical treatment is performed with high- 
intensity ultrasound waves. High-intensity ultrasound waves, when focused at a focal 
zone a distance away from the ultrasound source, will substantially medically affect 
tissue in the focal zone. However, the high-intensity ultrasound will not substantially 
affect patient tissue outside the focal zone, such as tissue located between the ultrasound 
source and the focal zone. Other treatment regimes of interest include unfocused high- 
intensity ultrasound wherein the ultrasound energy is distributed over a relatively broad 
region of tissue rather than being generally concentrated within a focal zone. 
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[0003] Ultrasound waves may be emitted and received by a transducer assembly. The 
transducer assembly may include a single transducer element, or an array of elements 
acting together, to image the anatomical tissue and to ultrasonically ablate identified 
tissue. Transducer elements may employ a concave shape or an acoustic lens to focus or 
otherwise direct ultrasound energy. Transducer such arrays may include planar, concave 
or convex elements to focus ultrasound energy. Further, array elements may be 
electronically or mechanically controlled to steer and focus the ultrasound waves emitted 
by the array to a focal zone to provide three-dimensional medical ultrasound treatment of 
anatomical tissue. In some treatments the transducer is placed on the surface of the tissue 
for imaging and/or treatment of areas within the tissue. In other treatments the transducer 
is surrounded with a balloon which is expanded to contact the surface of the tissue by 
filling the balloon with a fluid such as a saline solution to provide acoustic coupling 
between the transducer and the tissue. 

[0004] Examples of ultrasound medical systems and methods include: deploying an 
end effector having an ultrasound transducer outside the body to break up kidney stones 
inside the body; endoscopically inserting an end effector having an ultrasound transducer 
into the rectum to medically destroy prostate cancer; laparoscopically inserting an end 
effector having an ultrasound transducer into the abdominal cavity to destroy a cancerous 
liver tumor; intravenously inserting a catheter end effector having an ultrasound 
transducer into a vein in the arm and moving the catheter to the heart to medically destroy 
diseased heart tissue; and interstitially inserting a needle end effector having an 
ultrasound transducer into the tongue to medically destroy tissue to reduce tongue volume 
as a treatment for snoring. Methods for guiding an end effector to the target tissue 



WO 2005/058138 



3 



PCT/US2004/041799 



include x-rays, Magnetic Resonance Images ("MRI") and images produced using the 
ultrasound transducer itself, 

[0005] Low-intensity ultrasound energy may be applied to unexposed subdermal 
anatomical tissue for the purpose of examining the tissue. Ultrasound pulses are emitted, 
and returning echos are measured to determine the characteristics of the unexposed 
subdermal tissue. Variations in tissue structure and tissue boundaries have varying 
acoustic impedances, resulting in variations in the strength of ultrasound echos. A 
common ultrasound imaging technique is known in the art as "B-Mode" wherein either a 
single ultrasound transducer is articulated or an array of ultrasound transducers is moved 
or electronically scanned to generate a two-dimensional image of an area of tissue. The 
generated image is comprised of a plurality of pixels, each pixel corresponding to a 
portion of the tissue area being examined. The varying strength of the echos is preferably 
translated to a proportional pixel brightness. A cathode ray tube, computer monitor or 
liquid crystal display can be used to display a two-dimensional pixellated image of the 
tissue area being examined. 

[0006] When high-intensity ultrasound energy is applied to anatomical tissue, 
significant beneficial physiological effects may result. For example, undesired 
anatomical tissue may be ablated by heating the tissue with high-intensity ultrasound 
energy. By focusing the ultrasound energy at one or more specific focusing zones within 
the tissue, thermal effects can be confined to a defined region that may be remote from 
the ultrasound transducer. The use of high-intensity focused ultrasound to ablate tissue 
presents many advantages, including: reduced patient trauma and pain; elimination of the 
need for some surgical incisions and stitches; reduced or obviated need for general 
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anesthesia; reduced exposure of internal tissue; reduced risk of infection and other 
complications; avoidance of damage to non-targeted tissue; lack of harmful cumulative 
effects from the ultrasound energy on the surrounding non-target tissue; reduced 
treatment costs; minimal blood loss; and the ability for ultrasound treatments to be 
performed at non-hospital sites and/or on an out-patient basis. 

[0007] Ultrasound treatment of anatomical tissue may involve the alternating use of 
both low-intensity imaging ultrasound and high-intensity treatment ultrasound. During 
such treatment, imaging is first performed to identify and locate the tissue to be treated. 
The identified tissue is then medically treated with high-intensity ultrasound energy, such 
as for the purpose of destroying the tissue. After a period of exposure to high-intensity 
ultrasound, a subsequent image of the tissue is generated using low-intensity ultrasound 
energy to determine the results of the ultrasound treatment and provide visual guidance to 
the user to aid in subsequent treatments. This process of applying low-energy ultrasound 
to assist in guiding the position and focal point of the transducer, followed by high- 
energy ultrasound to ablate the undesired anatomical tissue, may continue until the 
undesired tissue has been completely ablated. 

[0008] In addition to imaging, monitoring of the temperature of the tissue being 
treated is desirable so that the tissue being treated, as well as the delivered treatment, can 
be readily visualized and controlled. For example, temperature monitoring is essential 
for hyperthermia treatments wherein tissue is exposed to ultrasound for the purpose of 
raising the temperature of the tissue. Hyperthermia treatments have been shown to be 
effective in the treatment of cancerous tumors, which are known to be more sensitive to 
heat than healthy tissue. Temperature monitoring is also useful for high-intensity 
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ultrasound treatment of tissue, such as ablation, since the temperature of the tissue 
provides an indication of the extent and spatial distribution of the treatment. 
[0009] The temperature of the tissue may be measured directly, such as with a 
thermocouple array. However, this arrangement is not desirable for several reasons. 
First, the thermocouple array must be placed into the region being treatment. This may 
require invasive placement of the array, thereby eliminating the non-invasive treatment 
advantages offered by ultrasound. In addition, the thermocouple array may have a 
relatively slow response time, providing the operator with delayed feedback regarding the 
status of the treatment. 

[0010] Various non-invasive temperature-monitoring methods have been developed to 
overcome the drawbacks of thermocouple arrays. For example, magnetic resonance 
imaging ("MRI") is used in the art to noninvasively measure temperature rises in vivo. 
However, MRI feedback requires an apparatus separate from the ultrasound treatment 
apparatus to provide the temperature measurement. As a result, MRI-assisted ultrasound 
methods are expensive and cumbersome, making such treatment impractical for routine 
use in surgical treatments. 

[0011] The prior art has suggested using ultrasound waves to measure temperature. 
Previous studies have shown that temperature rises in tissue cause local thermal 
expansion and sound speed changes, which can be estimated from perturbations of 
ultrasonic echo signals. Since the travel time of an ultrasonic echo depends on both the 
path length (which is changed by any thermal expansion) and sound speed (which 
depends on temperature), temperature changes in tissue being treated cause measurable 
alterations to echo travel times. By estimation and further processing of these travel 
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times, temperature rises can be approximately mapped. This arrangement is highly 
desirable, since heating effects of ultrasound treatment may be measured using 
inexpensive, portable, and unobtrusive ultrasound apparatus. For configurations in which 
therapy and monitoring are performed using the same ultrasound probe, this approach is 
particularly desirable, since a temperature map may be automatically co-registered with 
the therapeutic ultrasound beams. 

[0012] Available approaches to pulse-echo temperature mapping have employed 
cross-correlation of ultrasound echo signals to estimate changes in travel time (i.e., delay) 
for each pixel of an image, before and after heating. These travel-time maps are then 
smoothed and differentiated by signal processing to obtain maps of the echo strain, which 
is assumed to be proportional to the temperature rise of the tissue being treated. 
Repetition of this process can follow incremental temperature rises associated with 
prolonged heating. These approaches suffer from several drawbacks, including low 
spatial resolution and high artifactual content. Since many time-delay estimation steps 
are needed to build up a temperature image, these approaches are also computationally 
inefficient. 

[0013] There is a need for a method of measuring the temperature rise of tissue being 
treated by ultrasound that is noninvasive. There is a further need for method of 
monitoring temperature rise that is not expensive and is not cumbersome. There is a still 
further need for a method of monitoring temperature measurement that has relatively high 
spatial resolution and relatively low artifactual content, and is computationally efficient. 
Summary of the Invention 
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[0014] The present invention overcomes the limitations of the known art by estimating 
the time-dependent delay between two sets of narrow-band ultrasound signals. The travel 
time of a first signal is measured prior to treating anatomical tissue, then the travel time 
of a second signal is measured after the tissue is treated. The change in travel time 
between the first and second signals is used to generate a complete two-dimensional 
delay map from several two-dimensional multiplication and fast Fourier transform 
operations, rather than creating a map from a large quantity of individual time delay 
estimation steps. As a result, temperature maps can be generated more quickly as 
compared to prior methods. In addition, temperature maps generated using the present 
invention have fewer artifacts. 
Brief Description of the Drawings 

[0015] Further features of the present invention will become apparent to those skilled 
in the art to which the present invention relates from reading the following specification 
with reference to the accompanying drawings, in which: 

Fig. 1 A illustrates ultrasound signals before and after an ultrasound medical 

treatment; 

Fig. IB is a graph of changes in travel time for ultrasound echo signals; 
Fig. 1C is a graph of ultrasound echo strain; 

Fig. 2 is a flow diagram providing an overview of an ultrasound treatment 
method according to an embodiment of the present invention; 

Fig. 3 illustrates the relative amplitude and timing of ultrasound image frames 
and ultrasound treatments according to an embodiment of the present invention; 
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Fig. 4 is a flow diagram of a method for mapping temperature rise using pulse- 
echo ultrasound according to an embodiment of the present invention; 

Fig. 5 depicts the relative amplitude and timing of ultrasound image frames 
and ultrasound treatments according to an alternate embodiment of the present invention; 

Fig. 6 is a flow diagram of a method for mapping temperature rise using pulse- 
echo ultrasound according to an alternate embodiment of the present invention; 

Fig. 7 shows the relative amplitude and timing of ultrasound image frames 
according to another alternate embodiment of the present invention; 

Fig. 8 depicts a flow diagram providing an overview of an ultrasound treatment 
method according to another embodiment of the present invention; and 

Fig. 9 illustrates the relative amplitude and timing of ultrasound image frames 
and ultrasound treatments of the method of Fig. 8. 
Detailed Description 

[0016] It is well-known that the travel time of an ultrasound wave varies with the path 
length. The path length may vary in anatomic tissue due to thermal expansion of the 
tissue as a result of therapeutic ultrasound treatment, such as hyperthermia and ablation. 
Further, the travel time of an ultrasound wave varies with the speed of the ultrasound 
wave, which is a function of path temperature. Thus, changes in ultrasound wave travel 
times may be used to estimate temperature rise. 

[0017] Figs. 1A-1C provide an overview of a method for mapping temperature rise 
according to an embodiment of the present invention. Fig. 1 A generally illustrates a first 
ultrasound signal 1 1 returning to an ultrasound transducer (not shown) before treatment 
of a region of tissue, while ultrasound signal 13 represents an ultrasound signal returning 
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to the ultrasound transducer after treatment of the tissue. The relative delay between 
ultrasound echos returning to an ultrasound transducer, measured before and after 
treatment, is estimated as a function of depth, as shown in the graph of Fig. IB. The 
slope of the depth-dependent delay represents echo strain, as illustrated in Fig. 1C. Echo 
strain is defined as the slope (time derivative) of the time-dependent (equivalent to depth- 
dependent) signal delay, as shown in Figs. 1 A-1C. The amount of echo strain is 
proportional to the temperature rise of the heated region, as can be seen by comparing 
Figs. 1A-1C. 

[0018] A flow diagram of the general arrangement of an ultrasound treatment method 
10 according to an embodiment of the present invention is shown in Fig. 2. The method 
begins at step 12 by positioning proximate the tissue to be medically treated a transducer 
capable of transmitting and receiving ultrasound signals. At step 14 a low-intensity 
imaging ultrasound signal, such as a B-Mode signal, is emitted from the transducer and 
the reflected signals are collected to form a pre-treatment image frame "Fi." For each 
frame, a number of A-lines of raw echo signal radio frequency ("RF") data are obtained; 
a line number corresponds to azimuthal position and signal time corresponds to depth. It 
is understood that the terms "image" and "imaging" comprise, without limitation, any 
means for creating a visually-perceivable image that may be displayed, for example, on a 
monitor, screen or display, or printed. The terms "image" and "imaging" may further 
comprise image information that may be used by electronic devices such as computers 
without being displayed in visually perceivable form. 

[0019] Medical treatment begins at step 16 by emitting an ultrasound signal from the 
transducer. The term "medical treatment" as used herein may comprise several types of 
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ultrasound treatment. For example, the medical treatment may be a relatively low- 
intensity ultrasound signal, causing a small temperature rise in the tissue (such as 1 °C or 
less) that does not alter the tissue. Low-intensity ultrasound medical treatment is relevant 
for targeting of tissue for higher-intensity medical treatment, wherein the low-intensity 
ultrasound signal creates a small temperature rise that may be imaged to ensure that the 
desired region of tissue is targeted by the ultrasound beam. After targeting is confirmed, 
the amplitude of the ultrasound signal may be raised to a higher intensity. The high- 
intensity ultrasound signal may be used to administer a hyperthermia treatment or to 
ablate the tissue. 

[0020] After a portion of the tissue has been medically treated, a post-treatment image 
frame "F 2 " is generated at step 18. Information gathered from the pre-treatment and post- 
treatment image frames is then used at step 20 to compute the change in the temperature 
of the tissue being treated. The results of the temperature computations may then be 
displayed at step 22. The display may be in any convenient form to aid the user in 
interpreting the temperature data. For example, the display may be in the form of textual 
data relating to the temperature of the treated tissue. The display may also include a 
graphical presentation wherein colors are used to represent temperatures and temperature 
gradients of the treated tissue. If treatment is complete at step 24, the method is ended at 
step 26. However, if the tissue requires additional treatment, the transducer may be re- 
positioned if desired at step 28 using the image of step 18 as a guide. The method then 
returns to step 16 to again medically treat the tissue. The post-treatment image frame 
from the first treatment may be used as a pre-treatment image for the second treatment, as 
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shown by arrow 15 of Fig. 2. Alternatively, a new pre-treatment image 14 may be 
obtained, as indicated by arrow 17 of Fig. 2. 

[0021] Fig. 3 illustrates the method of Fig. 2 in relation to a time scale t. The B-Mode 
image scans 14, 18 are represented as frames Fi, F 2 respectively. For each frame, a 
number of A-lines of raw echo signal data are obtained, the number of each line 
corresponding to azimuthal position while signal time corresponds to depth. 
[0022] Referring now to Fig. 4 in combination with Fig. 3, a method for estimating 
the temperature rise of anatomical tissue in response to medical treatment according to an 
embodiment of the present invention is depicted. An ultrasound transducer is positioned 
proximate the anatomical tissue, then at step 14 a first image frame Fi (illustrated in Fig. 
4) is acquired. The frame may be stored electronically, such as in a computer, magnetic 
media and solid-state memory. The tissue may then be medically treated such as by 
hyperthermia or ablating it with high-intensity ultrasound waves, as at step 16. A second 
image frame F2 is then acquired at step 18 and electronically stored. For each frame, a 
number of A-lines of raw echo signal data are obtained. A line number assigned to each 
A-line corresponds to azimuthal position, whereas signal time corresponds to depth. The 
raw echo signals of frames Fj and F 2 may be processed at step 30, such as to obtain 
complex analytic signals by means of a conventional Hilbert transformation. For small 
temperature rises the echo signal p\ obtained after treatment differs frompo, the echo 
signal obtained prior to treatment, only by a time-dependent delay. Echo signals p 0 and 
pi may be expressed below in Equations 1 and 2: 

Po CO = «o (0 Pi W = a o (' + *) e~ icToit+so 

Equation 1 Equation 2 
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where <zq is the complex amplitude envelope, 0)q is a nominal center frequency of the 
signals, and St is the time-dependent delay. 

[0023] At step 32 the delay St is estimated as a function of time as follows. The 
processed echo signal p 0 from frame Fi is multiplied by the conjugate of p\ from frame 
F2. This conjugate product is spatially filtered to yield a spatial map of a zero-lag cross- 
correlation between Fi and F 2 . For the signals of Equations 1 and 2, a windowed zero-lag 
cross-correlation may be expressed by Equation 3 for a window centered at time t 0 : 



Since this zero-lag cross-correlation takes the form of a convolution between the 
conjugate product of the echoes and a windowing function w, this cross-correlation can 
be computed efficiently for all times to, using the fast Fourier transform. This eliminates 
the need for computing separate cross-correlations for multiple windows spanning the 
region of interest. 

[0024] The depth-dependent delay is then given by the phase of the zero-lag cross- 
correlation i? 01 (0,£ 0 ): 



The phase of the conjugate product is computed for all the A-lines, using techniques such 
as appropriate phase unwrapping, to efficiently estimate the azimuth- and depth- 
dependent delay St between the echo signals of the two frames using Equation 4. At step 
34 the relative strain between the two frames is computed as the derivative of this signal 




Equation 3 



arg{fl 01 (0,f)) 
co 0 



Equation 4 
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delay, taken along the time/range direction. This relative strain (termed an "echo strain 
map" herein) can be computed by known methods, such as finite differences. The echo 
strain map represents the slope of the depth-dependent delay and assumes a linear 
relationship between echo strain and temperature rise. 

[0025] At step 36, the echo strain map can optionally be spatially filtered to reduce 
any strain artifacts. Within a measurable constant of proportionality, the echo strain map 
can then be used at step 38 to estimate the spatially-dependent temperature rise between 
frames Fi and F2. 

[0026] The present invention provides several advantages over prior methods for 
strain estimation used in pulse-echo ultrasound temperature measurements. For example, 
the present invention is more efficient due to its relative simplicity, since single-filter 
operations may be performed on entire two-dimensional frames at once, rather than on 
numerous individually windowed signals as is the current practice. As a result, the 
present invention may be easily implemented for real-time measurement of anatomical 
tissue, using, for example, standard digital signal processing ("DSP") circuitry. In 
addition, the present invention provides temperature information having fewer artifacts 
and improved accuracy as compared to present methods. 

[0027] In a first alternate embodiment, multiple frames may be grouped together as a 
"frame set." A diagram of the relative amplitude and timing of ultrasound image frame 
sets and ultrasound treatments is depicted in Fig. 5. Frame sets FSi and FS 2 may be 
averaged prior to computing the conjugate product. In this embodiment the method of 
Fig. 4 may be utilized with minor modification wherein FS 1 , the average of the signals of 
the frames comprising frame set FSi, is substituted for Fi of step 14. Likewise, FS 2 , the 
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average of the signals of the frames comprising frame set FS2, is substituted for F 2 of step 
18. The resulting estimated temperature rise of step 38 contains a reduced number of 
uncertainties and artifacts as compared to the method of Fig. 4 using single frames. 
[0028] In a second alternate embodiment, the present invention may include averaging 
of echo conjugate products from multiple frames before computing and differentiating 
the phase map to obtain the echo strain map. A method according to this embodiment is 
illustrated by Fig. 5 in combination with Fig. 6. An ultrasound transducer is positioned 
proximate the anatomical tissue, then at step 14 a first image frame set FSi is acquired. 
The tissue may then be medically treated such as by hyperthermia or ablating it with 
high-intensity ultrasound waves, as at step 16. A second image frame set FS 2 is then 
acquired at step 18. The raw echo signals of frame sets FSi and FS 2 may be processed at 
step 30, such as to obtain complex analytic signals by means of a Hilbert transformation, 
as previously disclosed. At step 31 an average of the conjugate products of step 30 is 
calculated. The average may be computed, for example, by computing and then 
averaging echo conjugate products of: 

Fi ofFSiand Fi ofFS 2 ; 

F 2 ofFSi andF 2 ofFS 2; 

Fi ofFSi and Fi ofFS 2; 
and so on. The delay is computed at step 32, then at step 34 an echo strain map is 
generated using the averaged conjugate products. The resulting echo strain map is then 
spatially filtered at step 36, the temperature rise is estimated at 38, and the temperature 
map is displayed at 22, all in the manner previously detailed. The averaged conjugate 
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products serve to reduce uncertainties and artifacts in the estimated temperature maps as 
compared to the method of Fig, 4 using single frames. 

[0029] In a third alternate embodiment, a correlation coefficient may be computed by 
multiplying one of the complex frames by a phase compensation function e~ w , where 6 
is the low-pass filtered phase of the conjugate product. The normalized sum of the phase- 
corrected conjugate product is then a correlation coefficient whose amplitude indicates 
the reliability of the temperature estimation. Windowed sums, which can be 
implemented using known fast Fourier transform operations, provide a spatial map of the 
correlation coefficient. High values (i.e., values approaching 1) correspond to a high 
confidence in the temperature estimate, while values significantly less than 1 correspond 
to lower confidence. Estimated temperature maps can then be formed using higher- 
confidence estimates preferentially. 

[0030] In a fourth alternate embodiment of the present invention, the echo conjugate 
product can be computed using a two-dimensional filtering operation, in which a two 
dimensional (i.e., azimuth/array and depth/time directions) window is used instead of the 
one-dimensional window described above for the third alternate embodiment. This 
method of computing the conjugate product will result in smoother strain maps, and thus 
may reduce artifacts in the temperature images. 

[0031] In a fifth alternate embodiment of the present invention, the assumed 
relationship between echo strain and temperature rise may be non-linear rather than the 
linear relationship currently assumed. Such a non-linear relationship may be derived 
from measurements of the non-linear relationships between temperature, sound speed, 
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and thermal expansion. The non-linear relationship may also be obtained empirically 
from calibration measurements in tissues of interest. 

[0032] In a sixth alternate embodiment of the present invention, smoothing of strain 
maps may be enhanced by performing a thresholding operation before low-pass filtering. 
For example strain values above a predetermined temperature range of interest may be 
clipped. In this manner, artifactual strain peaks may be removed more effectively than by 
low-pass filtering alone. 

[0033] With reference to Figs. 7 and 8, in a seventh alternate embodiment of the 
present invention, the conjugate products of adjacent frames of a frame set FS may be 
obtained, computed and averaged. At step 16 the tissue is medically treated such as by 
hyperthermia or ablating it with high-intensity ultrasound waves. An image frame set FS 
is then acquired at step 18. The raw echo signals of frames Fi - F n may be processed at 
step 30, such as to obtain complex analytic signals by means of a Hilbert transformation, 
as previously disclosed. At step 31 an average of the conjugate products of step 30 is 
calculated. The average may be computed, for example, by computing and then 
averaging echo conjugate products of adjacent frames, i.e.: 

Fi andF 2 of FS; 

F 2 andF 3 of FS; 

F 3 andF 4 of FS; 

and so on. The delay is computed at step 32, then at step 34 an echo strain map is 
generated using the averaged conjugate products. The resulting echo strain map has 
reduced noise and artifactual content as compared to other imaging methods. The echo 
strain map is spatially filtered at step 36, the temperature change is estimated at 38, and 
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the temperature map is displayed at 22, all in the manner previously detailed. The 
temperature change occurring in the time interval separating adjacent time frames may 
thus be estimated and displayed. 

[0034] In an eighth alternate embodiment of the present invention the steps of the 
method of Fig. 2 may be generally followed for a series of medical treatments 16, shown 
as 16A-16C in Fig. 9. A pre-treatment image frame Fx is obtained at step 14 prior to the 
first treatment 16 A. A post-treatment frame F 2 is then obtained at step 18 and used with 
the pre-treatment image Fi to compute at step 20 the temperature change due to treatment 
16 A. A temperature map is displayed at step 22. Frame F 2 also serves as a pre-treatment 
image 14 for a subsequent treatment 16B, as shown in Fig 9. Likewise, frame F 3 is a 
post-treatment frame for treatment 16B and is also a pre-treatment frame for treatment 
16C. Finally, frame F 4 is a post-treatment image for treatment 16C. This process of 
imaging and treatment may be repeated as many times as desired. Further, the 
temperature change at 20 may be computed using any conventional mathematical 
function as previously discussed, such as averaging of frame images and averaging of the 
computed temperature change. 

[0035] The foregoing description of several expressions of embodiments and methods 
of the invention has been presented for purposes of illustration. It is not intended to be 
exhaustive or to limit the present invention to the precise forms and procedures disclosed, 
and obviously many modifications and variations are possible in light of the above 
teaching. It is intended that the scope the invention be defined by the claims appended 
hereto. 
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WHAT IS CLAIMED IS : 

1. A method for mapping temperature rise of anatomical tissue using pulse-echo 
ultrasound, comprising the steps of: 

a) obtaining a first signal of a first imaging ultrasound wave which has been 
reflected back from a region in the anatomical tissue at a first time; 

b) obtaining a second signal of a second imaging ultrasound wave which has 
been reflected back from the region in the anatomical tissue at a later second time 
wherein the tissue has received at least some medical treatment by the second time; 

c) computing first and second complex analytic signals from the first and second 
imaging signals; 

d) computing the depth-dependent delay from the conjugate product of the first 
and second analytic signals; 

e) generating an echo strain map from the slope of the depth-dependent delay; 

f) using the echo strain map to estimate the amount of temperature rise from the 
first imaging signal to the second imaging signal; and 

g) creating an image showing where temperature rise is occurring in the 
anatomical tissue. 

2. The method of claim 1, further comprising the step of spatially filtering the 
difference signal. 
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3. The method of claim 2, further including the step of performing a thresholding 
operation to remove artif actual strain peaks. 

4. The method of claim 1, further including the step of multiplying a signal of an 
ultrasound imaging wave by a phase compensation function to compute a correlation 
coefficient indicating the reliability of the temperature estimation. 

5. The method of claim 4 wherein the signal is the first imaging signal. 

6. The method of claim 5 wherein the signal is the second imaging signal. 

7. The method of claim 4, further including the step of implementing windowed 
sums to provide a spatial map of the correlation coefficient. 

8. The method of claim 7 wherein the windowed sums are two-dimensional. 

9. The method of claim 1, wherein the echo strain map representing the slope of 
the depth-dependent delay assumes a linear relationship between echo strain and 
temperature rise. 

10. The method of claim 1, wherein the echo strain map representing the depth- 
dependent delay assumes a non-linear relationship between echo strain and temperature 
rise. 
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1 1 . The method of claim 10 wherein the non-linear relationship is derived from 
measurements of non-linear relationships between temperature, sound speed and thermal 
expansion, 

12. The method of claim 10 wherein the non-linear relationship is derived 
empirically from calibration measurements in the anatomical tissue. 

13. The method of claim 1, wherein the medical treatment is ultrasound medical 
treatment. 

14. The method of claim 1, also including steps a) through g) for different 
regions to image the anatomical tissue, wherein the image includes medically-treated and 
medically-untreated regions of the anatomical tissue. 

15. A method for mapping temperature rise of anatomical tissue using pulse-echo 
ultrasound, comprising the steps of: 

a) obtaining a first set of frames comprising a plurality of imaging ultrasound 
wave signals which have been reflected back from a region in the anatomical tissue 
during a first period of time; 

b) obtaining a second set of frames comprising a plurality of imaging ultrasound 
wave signals which have been reflected back from a region in the anatomical tissue at a 
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later second time wherein the tissue has received at least some medical treatment by the 
second time; 

c) averaging together the signals of the first set of frames to obtain an averaged 
first imaging signal; 

d) averaging together the signals of the second set of frames to obtain an 
averaged second imaging signal; 

e) computing first and second complex analytic signals from the first and second 
averaged imaging signals; 

f) computing the depth-dependent delay from the conjugate product of the first 
and second analytic signals; 

g) generating a strain map from the slope of the depth-dependent delay; 

' h) using the echo strain map to estimate the amount of temperature rise from the 
first averaged imaging signal to the second averaged imaging signal; and 

i) creating an image showing where temperature rise is occurring in the 
anatomical tissue. 

16. The method of claim 15, further comprising the step of spatially filtering the 
difference signal. 

17. The method of claim 15, wherein the medical treatment is ultrasound medical 
treatment. 
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18. The method of claim 15, also including steps a) through i) for different 
regions to image the anatomical tissue, wherein the image includes medically-treated and 
medically-untreated regions of the anatomical tissue. 

19. A method for mapping temperature rise of anatomical tissue using pulse-echo 
ultrasound, comprising the steps of: 

a) obtaining a first set of frames comprising a plurality of imaging ultrasound 
wave signals which have been reflected back from a region in the anatomical tissue 
during a first period of time; 

b) obtaining a second set of frames comprising a plurality of imaging ultrasound 
wave signals which have been reflected back from a region in the anatomical tissue at a 
later second time wherein the tissue has received at least some medical treatment by the 
second time; 

c) computing complex analytic signals from a selected frame from the first set of 
frames and a selected frame from the second set of frames; 

d) computing the conjugate product of the complex analytic signals of step c); 

e) repeating steps c) and d) until conjugate products have been computed for all 
of the frames of the first and second frame sets; 

f) computing the average of the conjugate products of step e); 

g) computing the depth-dependent delay from the averaged conjugate product of 

step f); 

h) generating an echo strain map from the slope of the depth-dependent delay; 
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i) using the echo strain map to estimate the amount of temperature rise from the 
first averaged imaging signal to the second averaged imaging signal; and 

j) creating an image showing where temperature rise is occurring in the 
anatomical tissue. 

20. A method for mapping temperature change in anatomical tissue using pulse- 
echo ultrasound, comprising the steps of: 

a) obtaining a set of frames comprising a plurality of imaging ultrasound wave 
signals which have been reflected back from a region in the anatomical tissue during a 
period of time; 

b) computing complex analytic signals from the set of frames; 

c) computing the conjugate product of a pair of adjacent frames of the set of 

frames; 

d) repeating step c) until the conjugate products of all adjacent frames have been 
computed; 

e) averaging the conjugate products of the adjacent frames of step d); 

f) computing a depth-dependent delay map from the average conjugate product; 

g) generating an echo strain map from the slope of the depth-dependent delay; 

h) using the echo strain map to estimate the amount of temperature change from 
the first frame to the second frame; and 

i) creating an image showing where temperature change is occurring in the 
anatomical tissue. 
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21. A method for mapping temperature rise of anatomical tissue using pulse-echo 
ultrasound, comprising the steps of: 

a) obtaining a first signal of a first imaging ultrasound wave which has been 
reflected back from a region in the anatomical tissue at a first time; 

b) obtaining a second signal of a second imaging ultrasound wave which has 
been reflected back from the region in the anatomical tissue at a later second time 
wherein the tissue has received at least some medical treatment by the second time; 

c) computing first and second complex analytic signals from the first and second 
imaging signals; 

d) computing the depth-dependent delay from the conjugate product of the first 
and second analytic signals; 

e) generating an echo strain map from the slope of the depth-dependent delay; 

f) using the echo strain map to estimate the amount of temperature rise from the 
first imaging signal to the second imaging signal; 

g) creating an image showing where temperature rise is occurring in the 
anatomical tissue; and 

h) repeating the method at least once by: 

i) re-defining the second imaging signal obtained at step b) as the first 
imaging signal of step a); 

ii) obtaining a new second imaging signal at step b); and 
ii) repeating steps c) through g). 
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